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Abstract 

We introduce an algorithm for numerical calculation of deriva- 
tives of the Jones polynomial. This method gives a new tool for 
determining topology of knotted closed loops in three dimensions 
using computers. 
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1 Introduction 



In chemistry and biology, knotted macromolecules such as DNA 
knots can be synthesized in experiments. |2j In statistical mechan- 
ics of macromolecules, topological classification of configurations of 
polymer chains is an interesting and unsolved problem. [0, |5[] 
Topological constraints could severely restrict available degrees of 
freedom in the configuration space, and would have important ef- 
fects on thermodynamic quantities such as the entropy of the sys- 
tem. This problem was first formulated by Delbriick in '60s. 
Since then, several numerical works have been done. []7|, @, 0, [TDj In 
these works unknotting probabilities of closed loops are evaluated in 
different systems such as the closed self- avoiding random walks 
and the rod-bead model [OD 



In all these numerical researches, the special values (A^(— 1) or 
Ai4:(2)) of the Alexander polynomial Axit) have been used as tools 
for classifying knots. In the algorithms, a closed loop is considered 
as trivial if the special value of the Alexander polynomial is found 
to be one. However, it is known that the Alexander polynomial does 
not distinguish some knots from the trivial knot []TT|, |T2[ (see also 
T^]). Consequently, the special values of the Alexander polynomial 



are not complete invariants for classification. 

Recently various new link invariants, such as the Jones, the HOM- 
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FLY, and the Kauffman polynomials are introduced. [{TJ, [T5], |TB|, [T? 



If we could apply these new invariants to computer calculations, we 
would have a systematic method of determining topology of closed 
loops in 3 dimensions. However, it is not easy to calculate directly 
any special value of the new link polynomials for long closed loops 
with large number of steps. If we calculate the link polynomials 
through the Markov trace, the number of processes in computer cal- 
culation grows exponentially with respect to the step number of 
the loops. Furthermore, in numerical calculation of link polynomi- 
als, we have to assign a proper number to the variable t so that the 
polynomials do not become numerically too large (overflow). Thus 
it seems that we could not directly calculate the Jones polynomial 
for large closed loops using computers. 

The purpose of this paper is to show that if we devise new di- 
agrammatic method using oriented link diagrams, then derivatives 
of link polynomials can be easily calculated using computers. This 
approach has the following two advantages. 

1. We can evaluate derivatives of link polynomials in a calculation 
time proportional to some power of (the length of loops). 

2. We can avoid the overflow which could occur in evaluation of 
polynomials. 

The outline of this paper consists of the following. In §2 we in- 



troduce a state model for the Jones polynomials using oriented link 
diagrams, and construct an algorithm for calculating the derivatives. 
In §3 we show by examples that the derivatives of the Jones polyno- 
mial have enough information to distinguish knots from the trivial 
knot. In §4 applying the algorithm we estimate unknotting proba- 
bility of closed self-avoiding random walks. 



2 Algorithm of calculating the derivatives 
2.1 Oriented state model 

Let us consider diagrammatic methods for calculation of link polyno- 
mials. We call the methods state models. We have two types of state 
models of the Jones polynomial, that of unoriented link diagrams, 
and that of oriented link diagrams. We call a state model using ori- 
ented (unoriented) link diagrams oriented (unoriented) state model. 
L.H. Kauffman introduced an unoriented diagrammatic method called 
bracket polynomial for the Jones polynomial. It seems, however. 



that the bracket polynomial is not suitable for calculation of deriva- 
tives of the Jones polynomial using computers. 

We formulate an oriented state model of the Jones polynomial. 
Any oriented link diagram can be decomposed into the tangle 



diagrams given in Fig. 2.1. 
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Fig. 2.1 

For the tangle diagrams we define diagram weights. The nonzero 
diagram weights are given by the following. 

Uu = t\U2i = i,u'^ = i,u^' = r', (1) 

1) }} = G(+; l)i = 1, G(+; 1)11 = G{+; 1)S = t', 
G(+; l)g = 1 - t\ G(+; 2)}; = G(+; 2)i = 

2) 12 ^ 2)21 ^ 2)11 = - t\ 

Gi-; = G(-; l)i = 1, G(-; l)g = G(-; l)g = 

G(-; l)l\ = 1 - r^ G{-- 2)11 = G{-- 2)1 = r', 

G{-- 2)f, = G{-- 2f^ = 1, G(-; 2)}^ = _ ^-2. (2) 

The other nonzero diagram weights are given by the following rela- 
tions. 



Uab = Uab, U'' = U^\ 

G(±;l) = G(±;3), G(±; 2) = G(±; 4). 



(3) 



We introduce matrix notation X = X^^E^ 



Ei where (E-),, = 



6aiScj- For example we have 



Uab = Uab = 



1 



(4) 
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G(+;l) = G(+;3) 




^2 



(5) 



V 
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We now consider state sum of the oriented state model. We can 
represent any given link L by a link diagram. We decompose the 
link diagram into pieces of the oriented tangle diagrams given in Fig. 
2.1. We put variables on the edges of the tangle diagrams (see Fig. 



The variables can take two values, 1 and 2. Let each one of the vari- 
ables take either 1 or 2, then we have a configuration (of variables) 
on the link diagram. We now assign the diagram weights in eqs. (|I]) 
and (0) to the oriented tangle diagrams and take multiplication of 
the diagram weights. Then we take summation over all the possible 
configurations on the link diagram. Let the symbol (piL) denote the 
sum. We call the sum (j){L) oriented state sum. For example, we give 
the oriented state sum for the knot 3i (see Fig. 2.2). 



cPiL) = GfA-- 1) • Gfj,{--?>) ■ Glf{--^ ■ Uad ■ Uac ■ U^' ■ U-^'. (6) 



Here we have assumed the Einstein notation for sums over repeated 
indices. Another quantity necessary for calculation of the link poly- 
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2.2). 



Fig. 2.2 



nomial is the writhing number. [JT^ The writhing number is defined 



as follows. The oriented crossings have two types and G(— )). 

Each of these crossing types has been labeled with +1 or -1. The 
writhe of a link is defined by the sum of the labels for all the cross- 
ings in the link diagram. We thus define the Jones polynomial aj{L) 
using (j){L) and writhe by 

aj{L) = . {t'yrUhe^ 



t^ + t 



By checking the Reidemeister moves of the tangle diagrams []TB 



we 



can show that the definition (|7|) gives an isotopy invariant of knots 
and links. It is the Jones polynomial. 

The oriented state model has the advantage that all the nonzero 
entries of the diagram weights of U and U do not involve the ± 
sign although the nonzero entries are off-diagonal. If we take the 
limit t ^ 1, then the oriented state model satisfies the following 
properties. 

limU = lim U = , 
t^i t^i \l oj 

\iu,G-\{k-±) =5^,51 (for A: = 1,2,3,4). (8) 

where (5^ and (5^ denote the Kronecker delta. 

Let us compare the oriented state model with the unoriented one. 
If we use the unoriented state model (bracket polynomial), in the 



limit t ^ 1, the nonzero entries of the matrices U have ± factors. 
Therefore we have to count how many times the signs ± occur in a 
given configuration, and thus calculation time will grow as fast as 2^ 
with respect to the system size A^. 



2.2 Expansion of the state sum 

Let us discuss calculation of derivatives of the Jones polynomial. We 
set 

t^ = l + e, (9) 

and consider expansion of the Jones polynomial with respect to e. 
Because the Jones polynomial is link invariant, the coefficients are 
also link invariants. The oriented state model in the last section 
directly leads to an algorithm for calculation of the coefficients. We 
expand the matrices G and U in terms of the parameter e. We 
calculate the oriented state sum using the expanded matrices. For 
an illustration we consider the matrices G{+; 1) and Uab 

^0 0^ 
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We call the matrices of 0-th, 1-st, 2-nd orders • • • in the e expansion 
matrices of 0-th order, 1-st order, 2-nd order, • • •, respectively. 

Let us show how to evaluate the second order contributions in the 
oriented state sum for the knot 3i (Fig 2.2). We first expand 0(i^) 
in eq. (g) by e. 

4>{L) = 2 



+ e 




+ Gfji-; l)msiSp{f{-, 4)(i)f/i3'f/(°'c7°(»)c7^-''(0) + . . .} 
+ o(e\ (11) 

We now consider the first term of the 2-nd order Ggj(— ; l)*^^)- 
Gg'l{—; 3)^^) • • • in the expansion (|TT]). This term is equivalent to the 
configuration of Fig 2.3 where the two matrices enclosed with circles 
are 1-st order and all the other are 0-th order. 

We can readily evaluate the term without directly calculating the 
oriented state sum with respect to all the variables on the diagram. 
We identify the two variables on those edges of the braiding diagrams 
that are connected by the dotted lines each other, and then take the 
sum for the remaining independent variables. For any link diagram 
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we can thus reduce the number of independent variables into at most 
4 in calculation of the second order contributions. Thus we have the 
following 

a,b,c,d,e,f,g,h 

= Y. 'G°}(-;l)(').G5l(-;3)<') = -4. (12) 

a,b,g,h 

Here the symbol Y^abgh means that we take the summation for the 
independent variables a, 6, h where c, (i, e and / are assumed to be 
conjugate of h^a^g and /i, respectively; when a = 1, then d = 2. It 
is important to note that in calculation of the term in (|T2|) we can 



eliminate the parts of the Kronecker deltas and the matrices U^^^ 
because of the property (||). 

To finish this section, we make a comment. We can apply the 
method in this section to other link polynomials. We can calculate 
derivatives of the HOMFLY polynomial, in particular, of the Alexan- 
der polynomial. For an illustration, the oriented state model for the 
Alexander polynomial is given in Appendix A. 

3 On derivatives of the Jones polynomial 
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3.1 Tools of classifying knots 

We give a table of the derivatives of the Jones polynomial and the 
special value of the Alexander polynomial (A/^ (— 1)) for knots. 

Table 3.1 

It is easy to see from the table that derivatives of the Jones polyno- 
mial are useful in classifying knot types. Since link polynomials are 
Laurent polynomials of a variable we can recover link polynomials 
from their derivatives (if we previously know something about the 
links such as upper bound of the crossing number etc.). 

You may think that the second derivative of the Jones polynomial 
is not sufficient to distinguish a non-trivial knot from different knots. 
However, if we calculate the third order, then we may have much 
more precise distinction of knots (at least as far as prime knots are 
concerned). Further, if we calculate the fourth order we can distin- 
guish even such knots that cannot be distinguished by the Alexander 
polynomial. 

Let us consider a composite knot (3i)(3i)(3i) and a prime knot 8io 
that share the same Alexander polynomial (see e.g., Ref. [[TOj)- From 
Table 3.1 we see that the fourth derivative of the Jones polynomial 
distinguishes the two knots while for the Alexander polynomial at 
t = —1 {Ak{—1)) the values for knots (3i)(3i)(3i) and 8io coincide. 
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When we calculate the Alexander and the Jones polynomials for 
composite knots (links) it is useful that the link polynomial for the 
composite link Li^L2 is given by a{Li^L2) = a{Li)a{L2). 

We think that the algorithm can make up for possible faults of 
the Alexander polynomial and thus we can determine knot types of 
closed loops much more exactly. 

3.2 Remarks on derivatives of the Jones polynomial 

It is noteworthy that derivatives of the Jones polynomial have many 
applications in knot theory. The following formulas have been given. 

jVK{t =1) = 0, (13) 

—rVKit = 1) = const. X a2(K). (14) 

Here a2{K) is the second coefficient of the Conway polynomial (the 
Alexander polynomial) [|r2j , which is related to Arf invariant. [jT^, |21[] 
We may discuss derivatives of the Jones polynomial from the view- 
point of perturbation theory of the Chern-Simons field theory. The 
coefficients for the terms in the perturbational expansion of the Wil- 
son lines are expressed in terms of the Feynman integrals. []22| The 



numerical calculation of derivatives of the Jones polynomial corre- 
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sponds to evaluation of the Feynman integrals. 

4 Numerical results and concluding remarks 

We can apply the algorithm introduced in the present paper to any 
system of closed loops. For an illustration using the second derivative 
of the Jones polynomial we analyzed the system of rings [jS] generated 
by closed Gaussian random walks. 

Graph 4.1 

The probability of occurrence of the trivial knot is plotted in the log- 
arithmic scale as a function of the number of steps (A^) of generated 
closed loops. We call the probability of the trivial knot unknotting 
probability. With the number of steps fixed we evaluate the un- 
knotting probability Pq = Nt/M by counting the number A^^ of trivial 
knots when we generate M= 1000 closed random walks. 

From Graph 4.1 we see that the unknotting probability has an 
exponentially decaying behavior Po{N) ~ exp(— A^/A/'q), where we 
call A^o characteristic length. We may consider A'o necessary steps to 
form a nontrivial knot. From Graph 4.1. we have an estimation of 
characteristic length A'o; 300 < A'o < 370. 

Let us consider efficiency of our algorithm. We plot the time 
necessary for calculation of the second coefficients as a function of 
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the number of steps (A^) of given closed loops. 

Graph 4.2 

We see from Graph 4.2 that the calculation time for the second 
derivative of the Jones polynomial behaves asymptotically as N'^ 
with respect to N. From the result we may expect that the time 
for calculation of the r-th derivative grows as N'^, or at least non- 
exponentially with respect to N. 

This result makes a clear difference from the other approaches 
to the Jones polynomial. We can calculate the Jones polynomial 
through explicit evaluation of the Markov trace or diagrammatic 
calculation using the skein relation. If these methods are applied to 
computer calculations, then they yield exponentially growing calcu- 
lation time with respect to A^. We thus conclude that derivatives of 
link polynomials give a simple and systematic computer-orientated 
method of determining topology of given closed loops in 3 dimen- 
sions and that the algorithm in the paper can be applied to various 
problems related to knotted configurations in statistical physics and 
many body problem. 
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A Appendix 



We introduce an oriented state model of the Alexander polynomial. 
TD| The diagram weights Uaf ■ •, and (^(zb; /c)"^ are defined for 



the tangle diagrams in Fig. 2.1. The nonzero diagram weights are 
given by the following. 

Uu = U'' = Un = U'' = U22 = U^^ = 1, U22 = = -1, (A.l) 

G(+; 1)}} = 1, G(+; l)g = G(+; 1)11 = G{+- 1)1 = 1 - r^ 
G(+; 1)1 = -t-\ G(+; 2)}} = 1, G(+; 2)g = G(+; 2)11 = 
G(+; 2)?? = - 1, G(+; 2)1 = -t-\ G(+; 3)}} = 1, 
3)g = 3)?^ = 3)}^ = 1 - t-\ 

3)1 = -t-\ 4)}} = 1, 4)g = 4)?! = r^, 
G(+; 4)11 ^ ^ _ ^-4^ 4)22 ^ _^-4^ ^ 1^ 

l)g = G{-- l)?i = l)g = 1 - 

1) 1 = -t\ G{-- 2)}} = 1, G(-; 2)g = G(-; 2)g = t\ 

2) 22 = — 1, G{—] 2)22 = — ^"^5 ; 3)1} = 1, 

3) 21 = G{—] 3)^2 ~ G{—] 3)2! = 1 — t^, ; 3)22 = — 
G(-; 4)}} = 1, G{-- 4)g = 4)^ = t\ G{-- 4)?? = I - t\ 
G{-Afi = -t\ (A.2) 



By the oriented state model we can calculate derivatives of the 
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Alexander polynomial. We give two remarks. (1) U"'^ • • •, are diago- 
nal matrices, although there are —1 in some of the nonzero entries. 
(2) The nonzero off-diagonal elements of the matrices G(± : j) be- 
come 1 when t ^ 1 (cf. eq. (§)). 
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Figure Captions 

Fig. 2.1 Oriented tangle diagrams. 

Fig. 2.2 Link diagram of knot 3i and variables a, 6, c, (i, e, /, /i, j. 
Fig. 2.3 Configuration for the term 

a,6,c,(i,e,/,5i,/i 

Table 3.1 The symbol Ak{—1) denotes the Alexander polynomial. 
Remark that the first derivative of the Jones polynomial for a knot 
vanishes (see the formula (O))- 

Graph 4.1 Unknotting probability Po{N) as a function of length A^. 
Graph 4.2 Calculation time versus length N. 
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